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The features of the IV characteristics of metallic carbon nanotubes (m-NTs) in different exper- 
imental setups are studied using semi-classical Boltzmann transport equation together with the 
heat dissipation equation to account for significant thermal effects at high electric bias. Our model 
predicts that the shape of the m-NT characteristics is basically controlled by heat removal mech- 
anisms. In particular we show that the onset of negative differential resistance in free standing 
nanotubes finds its origins in strong transport nonlinearities associated with poor heat removal 
unlike in substrate-supported nanotubes. 
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Since their discovery Q] , carbon nanotubes (NTs) have 
captured the attention of both the scientific and tech- 
nological communities because of their mechanical sta- 
bility as well as their high thermal and electrical capa- 
bilities Q not usually seen in other materials. Among 
their outstanding properties is the peculiarity to behave 
as metals or semiconductors depending on their chirality 
0, Q ■ Their capability of carrying large current densi- 
ties at room temperature make them prominent mate- 
rials for field effect transistors [fj and interconnects 
in high speed nanoscale electronics. The first mea- 
surements of the IV characteristics on individual metal- 
lic nanotubes (m-NTs) were performed on substrate- 
supported tubes in a configuration similar to a field effect 
transistor. The output characteristics resembled those 
obtained in semiconductors as they exhibited linear de- 
pendence on voltage in the low field regime, and satura- 
tion at 7~25/iA under high biases. Later measurements 
have shown that in short m-NTs, the saturation level 
could be overcome, but a nonlinear behavior still persists 
in substrate-supported configurations 0, 0] . In recent 
experiments performed on m-NTs standing freely across 
a trench, current levels in the high bias regime were a 
few times smaller than those in nanotubes supported by 
substrates while negative differential resistance (NDR) 
was reported [Tl| . This behavior, indicating significant 
heat production/dissipation in m-NTs devices, was inter- 
preted in terms of the onset of a non-equilibrium optical 
phonon population. 

In this Letter we show that the nonlinearities in the 
IV characteristics of m-NTs, find their origins in Joule 
heating and the efficiency in heat removal from the HI- 
NT. We determine the temperature profile along nan- 
otubes by solving the Boltzmann transport equation si- 
multaneously with the heat transfer equation. We specif- 
ically demonstrate that while the nonlinear behavior in 
substrate-supported tubes (SSTs) (Fig.^a) emerges due 
to a strong imbalance between carrier distribution func- 



tions with positive and negative Fermi velocity, the char- 
acteristics in free-standing tubes (FSTs) (Fig.Qb) essen- 
tially arise from an inhomogeneous self- heating effect. 

The distance between sub-bands in m-NTs using tight- 
binding calculations can be estimated as AE met ~ 
67oa/D where a = 0.14nm, 70 ~ 3eV and D is the di- 
ameter of the tube Hence, for transport purposes, 
the electronic structure of small diameter m-NTs close to 
the Fermi level can be well described by linear dispersion 
relations E{k) = ±hvp(k — kp) with a Fermi velocity vf 
close to 8.10 7 cm/s. This approximation is valid as long 
as temperature satisfies that fc^T <^AE met and the bias 
voltage V < AE met . The current along the nanotube is 
naturally defined as: 
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where n + (n 



is the electron density with positive (neg- 
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ative) Fermi velocity. Similarly to Ref. |13| . we assume 
that carrier populations are described by Fermi statistics 
with different quasi- Fermi level E F depending on the sign 
a of the Fermi velocity. Using the zeroth moment of the 
Boltzmann equation, the electric field F satisfies: 

v F d x (n++n-)-^ = 2C ph (I, T). (2) 

Here C p h{l , T) is the momentum integral of the collision 
integral considering only the contribution to scattering 
involving high energy optical phonons (huj ~ 0.18eV). In 
our model, phonons are assumed to be in local thermal 
equilibrium with electrons, and their occupation number 
follows the Bose-Einstein statistics. 

A realistic description of the IV characteristics re- 
quires the consideration of the temperature profile T(x) 
along the m-NT in the right hand side of Eq. [21 We ac- 
count for this local temperature variation by solving the 
heat production/dissipation equation: 



d / dT\ + 
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in which k, 7, T , and q* are the thermal conductiv- 
ity, the thermal coupling to the substrate, the substrate 
temperature, and the heat dissipation per unit volume, 
respectively. The first term on the left hand side of Eq. |3| 
corresponds to the heat flowing through the leads; the 
second, to the flow driven through the substrate or any 
surrounding medium. The heat produced locally is as- 
sumed to obey Joule's law q* = IF/irDt, being t the 
effective thickness of the tube (t rj 0.34nm). The HI- 
NT thermal conductivity dependence on temperature has 
been predicted to follow k(T) = krtTrt/T , due to the 
Umklapp process of acoustic phonons [fj] for T>"iQQK. 
In this equation, krt is the thermal conductivity at a 
reference temperature Trt i K ~ 20W/cmK at 300K). 
An estimate of the thermal coupling constant 7 of nan- 
otubes standing on SiC>2 substrates is 10 1:L W/cm 3 K [l3j . 
The explicit dependence of the temperature on the local 
field requires Eqs. [21 and [3] to be solved simultaneously, 
leading to a non-linear second order differential equation 
for the temperature. We set the boundary conditions as 
T(±L/2) = To, where L is the m-NT length, assuming 
that both contact leads have the same temperature as 
the substrate (i.e. we neglect the heating at the leads 
due to a contact resistance). 
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FIG. 1: Experimental setups: (a) Substrate-supported nan- 
otube (SST); (b) Free-standing nanotube (FST). (c) IV char- 
acteristics for 800nm-long nanotube with different values of 
coupling coefficient 7. (d) Temperature in the middle of the 
tube AT(0) as a function of the current I for different 7. 

In Fig.^c) we compare the performances of an 800nm 
long m-NT assuming different 7 values. The parameters 
used in the calculations are k = 20W/chlK, t = 22fs and 
hid = O.I8eV. For the smallest 7 values, the m-NT ther- 
mal coupling to the substrate is vanishingly small (less 
than 10% of the power generated along the tube is dissi- 
pated in the high bias through the substrate) and NDR 
is observed in the IV characteristics. As we increase 7 to 
~ 10 11 W/cm 3 K saturation is observed. If we further in- 
crease 7, the saturation is overcome but deviations from 



the linear regime still arise. Our results clearly indicate 
that in the low bias regime, thermal effects can be ne- 
glected, but the high bias regime is basically governed by 
the m-NT capability for heat removal. Hence, the fea- 
tures of the electrical characteristics depend not only on 
the tube's length but also on 7. In Fig. ^d) we plot the 
temperature difference AT(0) = T(0) - T at the m-NT 
midlength as a function of the current for the set same of 
7- values. Our model predicts that the larger the coupling 
coefficient 7, the higher current in the m-NT. 
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FIG. 2: Features comparison of 800nm long m-NT strongly 
(left) and weakly (right) coupled to the substrate: (a) IV 
characteristics; (b) Temperature profile AT(x) as a function 
of the voltage V; (c) Fermi distribution at midlength for each 
of the carrier branches at different points indicated in the IV 
characteristics. The distance between dashed line corresponds 
to the optical phonon energy. 

In Fig. [21 we consider the situation of an 800nm 
long m-NT in two thermal coupling regimes i.e. i) 
strongly coupled (SST with 7 = 10 12 W/cm 3 K, left hand 
side of Fig. [21 and ii) weakly coupled (FST with 7 = 
10 9 W/cm K, right hand side of Fig. [21 to the substrate. 
In the SST, the TV characteristic exhibits a nonlinear be- 
havior with high resistance at high fields (Fig. 0a. LHS). 
In this case, the substrate acts as a heat sink, provid- 
ing efficient dissipation of the heat produced along the 
tube while keeping both the electric field and power dis- 
sipated per unit volume low. This efficient heat removal 
reduces the temperature rise in the m-NT, thereby fa- 
voring relatively high current levels while maintaining a 
quasi uniform temperature profile along the tube (ex- 
cept at the edges) as shown in the LHS of Fig. [2]b. The 
LHS of Fig. [2Jc displays the profile of the (Fermi-like) 
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carrier distributions at the NT mid-length for electrons 
with positive and negative Fermi velocity (/ + and / _ ) 
at the different points indicated along the IV curve in 
Fig. [2|a (LHS). No significant heating is observed before 
the separation between the two distribution quasi-Fermi 
levels becomes comparable to the optical phonon energy 
(Ep~Ep Kihuj). Hence, the nonlinear behavior observed 
in this kind of configuration is basically attributed to the 
smooth onset of optical phonon scattering || [Jj| , despite 
minor thermal broadening induced by Joule heating on 
the carrier distributions. In the FST, the IV characteris- 
tic exhibits a NDR at relatively low biases (Fig. [21a RHS), 
and the current level is lower than in the previous case of 
strong thermal coupling. Because of poor heat removal, 
temperature rises rapidly with external biases manifest- 
ing a pronounced maximum at NT mid-length (RHS of 
Fig. [21b) . One notices the different scales on the volt- 
age and temperature axis between the SST (LHS) and 
FST (RHS). Simultaneously, the Fermi distributions /+ 
and / _ remain weakly separated but experiencing signif- 
icant heating even at low voltages as shown in Fig. [2Jc 
(RHS) for the set of points along the IV characteristic 
on Fig. [2a (RHS). The fast rise of temperature, conse- 
quence of the poor heat removal, locally broadens the 
electron distributions, and therefore enhances the scat- 
tering locally even at low bias. The rapid onset of non- 
uniform T-profile induces nonlinearities in the electric 
field, which result in a current diminution in order to 
satisfy the boundary conditions on the applied voltage 
at the NT contacts. 
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FIG. 3: Comparison between the IV characteristics for FSTs 
of different length calculated with our model (lines) and ex- 
perimental data (symbols) of Ref. [Til ]. 

In Fig. [3 we compare the results of our model for FSTs 
with the experimental data of Ref. , for a which good 
agreement is observed for all NT-lengths. The physi- 
cal parameters used in the calculations are ?ia; = 0.I8cV 
and krt = 20W/cmK. The relaxation time (r ~ 30fs) 
used to compute the characteristics increases with diame- 
ter consistently with theoretical predictions 0] , and lies 



within previous experimental estimates |1 0, [Jj| ■ Our 
model also assumes the presence of a contact resistance 
{R c ~ 20kf2). 

In conclusion our model shows that the high field trans- 
port properties of m-NTs are strongly controlled by the 
onset of thermal effects, which can be altered by modify- 
ing the experimental setup. Using realistic parameters, 
our model is able to reproduce the IV characteristics in 
both substrate-supported and free-standing nanotubes. 
However the nonlinear behavior observed in each of the 
systems is caused by different phenomena affecting the 
carrier distribution. In the former, the transport features 
are due to an imbalance in the carrier population with 
positive and negative Fermi velocities associated with 
strong scattering and high current levels. In the latter, 
the NDR is attributed to strong scattering enhancement 
due to the broadening of the carrier distribution with low 
current level. 
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